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Abstract 

We apply the theory of random Schrodinger operators to the analysis of multi-users commu- 
nication channels similar to the Wyner model, that are characterized by short-range intra-cell 
broadcasting. With H the channel transfer matrix, HH^ is a narrow-band matrix and in many 
aspects is similar to a random Schrodinger operator. We relate the per-cell sum-rate capacity 
of the channel to the integrated density of states of a random Schrodinger operator; the latter 
is related to the top Lyapunov exponent of a random sequence of matrices via a version of the 
Thouless formula. Unlike related results in classical random matrix theory, limiting results do 
depend on the underlying fading distributions. We also derive several bounds on the limiting 
per-cell sum-rate capacity, some based on the theory of random Schrodinger operators, and some 
derived from information theoretical considerations. Finally, we get explicit results in the high- 
SNR regime for some particular cases. 
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I. Introduction 

The growing demand for ubiquitous access to high-data rate services, has produced a 
huge amount of research analyzing the performance of wireless communications systems. 
Techniques for providing better service and coverage in cellular mobile communications 
are currently being investigated by industry and academia. In particular, the use of joint 
multi-cell processing (MCP), which allows the base-stations (BSts) to jointly process their 
signals, equivalently creating a distributed antenna array, has been identified as a key tool 
for enhancing system performance (see [1][2] and references therein for surveys of recent 
results on multi-cell processing). 

Motivated by the fact that mobile users in a cellular system "see" only a small number 
of BSts, and by the desire to provide analytical results, an attractive analytically tractable 
model for a multi-cell system was suggested by Wyner in [3] (see also [4] for an earlier 
relevant work). In this model, the system's cells are ordered in either an infinite linear 
array, or in the familiar two-dimensional hexagonal pattern (also infinite). It is assumed that 
only adjacent-cell interference is present and characterized by a single parameter, a scaling 
factor a G [0,1]. Considering non-fading channels and a "wideband" (WB) transmission 
scheme, where all bandwidth is available for coding (as opposed to random spreading), 
the throughputs obtained with optimum and linear MMSE joint processing of the received 
signals from all cell-sites are derived in [3]. Since it was first presented, "Wyner-like" models 
have provided a framework for many works analyzing various transmission schemes in both 
the up-link and down-link channels (see [1] [5] and references therein) . 

In this paper we consider a generalized "Wyner-like" cellular setup and study its per-cell 
sum-rate capacity. According to Wyner's setup, the cells are arranged on a circle (or a line), 
and the mobile users "see" only a fixed number of BSts which are located close to their cell's 
boundaries. All the BSts are assumed to be connected through an ideal back-haul network 
to a central multi-cell processor (MCP), that can jointly process the up-link received signals 
of all cell-sites, as well as pre-process the signals to be transmitted by all cell-sites in the 
down-link channel. The model is characterized by short-range intra-cell broadcasting. Thus, 
if we denote by H the channel transfer matrix, then HH^ is in many aspects similar to a 
random Schrodinger operator. More specifically, the per-cell sum-rate capacity of the channel 
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is a function of the integrated density of state of HH\ which in turn is related to the top 
Lyapunov exponent of a random sequence of matrices via a version of the Thouless formula. 
Unlike associated results in classical random matrix theory, limiting results do depend on 
the underlying fading distributions. 

As an application of our result and motivated by the fact that future cellular systems 
implicitly assume high-SNR configurations mandatory for high data rate services, we get 
explicit results in the high-SNR regime for some particular cases. 

The rest of the paper is organized as follows. In Section [Til we present the problem 
statement. In Section UTTl we prove the convergence of the per-cell sum-rate capacity when 
the number of cells and BSts goes to infinity and we express the limit in terms of the 
Lyapunov exponent of a sequence of random matrices (Theorem |2]). In Section ITVl we give 
several reformulations of this result that yields a particularly simple expression in the high- 
SNR regime. In Section |V] we give different bounds on the per-cell sum-rate capacity, some 
of which are based on the theory of product of random matrices, and some on information 
theoretical considerations. In Section IVI| we specialize the results and make them explicit in 
some particular cases. Finally in Section IVHI we discuss some open problems using numerical 
simulations. The relevant background on the theory of Lyapunov exponents is given in 
Appendix IA.ll and the relevant background on exterior products is given in Appendix ICl 
Several proofs are postponed to Appendices IA.21 IA.3I and El The per-cell sum-rate capacity 
of the non-fading channels is derived in Appendix [Dl 

II. Problem statement 

In this paper we consider the following setup, m + d cells with K single antenna users per 
cell are arranged on a line, where the m single antenna BSts are located in the cells. Starting 
with the WB transmission scheme where all bandwidth is devoted for coding and all K users 
are transmitting simultaneously each with average power p, and assuming synchronized 
communication, a vector baseband representation of the signals received at the system's 
BSts is given for an arbitrary time index i by 

y(i) = H m (i)x(i) +n(i), 

where x(i) is the (m + d)K complex Gaussian symbols vector, z(i) is the unitary complex 
Gaussian additive noise vector. Note that the SNR is p. From now on, we omit the time 
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index i. H m is the following m x Kim + d) channel transfer matrix, which is a d + 1 block 
diagonal matrix defined by 

Cm+i o ••• \ 
C24+1 C24+2 '■ 



Cm,m Cm,m+1 ' ' ' Cm,d+m J 

where Qj are 1 X K row vectors. For s £ N*, we will denote by ( s the vector (Ca-d,a> • • • j Cs,«) 
and we denote by 7r it distribution. We assume in the rest of the paper that for n £ N* and 
< % < d the vectors (C n -i,n) are distributed according to 7Tj. We define Q = (C n )neN* and 
P, the probability distribution on Q associated to the above problem. We denote by E the 
associated expectation. We also use the 2 norm for vectors and matrices. For matrices, it is 
the Froebenius norm, which is a sub-multiplicative norm. 

Throughout this paper, we assume a subset of the following hypotheses. 
(HI) The vectors (OOjeN* form a stationary ergodic sequence. 
(H2) There exists e > such that for < i < d, K ni |log |a;|| 1+e < 00. 
(H3) If (x , . . . ,Xd) is distributed according to ir, then almost surely, x x\ ^ 0. 

For m £ N* and A > 0, we set G m = HmH^ + AId m , where Id m is the m x m identity 
matrix. Although G m depends on A, we will not write that dependence unless there is an 
ambiguity. Under the assumption that H m (i) is ergodic with respect to the time index i, that 
the Channel State Information (CSI) is known at the receiver whereas the users know the 
statistics of the CSI, and that the channel varies fast enough so as to allow each transmitted 
codeword to experience a large number of fading states, we follow [1] and study the per-cell 
sum-rate capacity that is given by the following formula ([6]) 

Ca Pm (p) = -Elogdet (ld+ P H m Hl) = logp + — E (log det G m (X)) , (1) 
m m 

where A = 1/p. 



(Cl,l Cl,2 
C 2 ,2 
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III. Main result 



We set for i e N* 



( 



a 



Cd(t-l)+l,d(i-l)+l Cd(t-l)+l,d(i-l)+2 
Cd(i-l)+2,d(i-l)+2 



(d(i-l)+l,d^ 
Cd(i-l)+2,di 



and 



D; 



V 



^ Cd(i-2)+l,d{i-l)+l Cd(i-2)+2,d{i-l)+l 
Cd(t-2)+2,d(i-l)+2 



di.di 



J 



Cd(t-l),d(t-l)+l 
Ci(t-l),£i(t-l)+2 



V 



o •■■ o Cd{i-i),di J 

For all i 6 N*, Cj are d x dK matrices and Di are dK x d matrices. We fix with z < 
or j < so that CiDi = Id^d- 

We thereby get the following block description of H 

\ 



H, 



dn 





Di 


Od,dK 




Od,dK 


Od,dK 


c 2 


Di 














Od,dK 


Od,dK 




Od,dK 


G n 


Di +l 



J 

where Od,dK is the d x dK zero matrix. 

Under the hypothesis (HT5J), in order to study the limit in m of Cap m (p), it is enough 
to study Cap n d{p) (see Remark [321 following the proof of Lemma [25]). We get the following 
block representation of Gdn- 



G 



dn 



( C-SA + D\D 2 + A Id d 


{C 2 D 2 ) ] 


o d 


o d 


C 2 D 2 






o d 








(c n D n y 


V o d 




C n D n 


C n Cl + Dl +1 D n+1 + XId d 



Note that under (HH]), for all i e N*, CiDi is a d x d invertible matrix. 
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For i G N*, we denote by Mj the following matrix 







-(Ci+iA+O^CiA 


-(Ci+iA+i) - ^ (acj + ^I+iA+i + AId d ) 



and denote iVj = /\ Mj. Moreover, 7(iV) denotes the top Lyapunov exponent associated 
with {iVi}, i.e. 

7(JV) = lim —log ||iV„ — iVill . 

n^oo Tl 



Note that by Theorem [20], 7(iV) is deterministic. See Appendix IA.1I for the definitions 
concerning the Lyapunov exponents and AppendixOfor the relevant background on exterior 
products. Recall that (Mj) ieN * and (iVj)j 6 N* depend on A. 



Theorem 2 Assume (Iflty, and (HB\), and set A = 1/p. 

1. VKe /jai>e 

Cap m {p) ► log p + E n log 



CoC]| + ^7 (AO =C>(p), 



where the expectation is taken such that (( , . . . , is distributed according to it. 
2. As p goes to infinity, 

Cap(p) = logp + E w log |CoC]| + \l (N{\ = 0)) + o(l). 

The theorem is proved in Appendix [A] 

As an alternative to deriving exact analytical results we will also be interested in extracting 
parameters that characterize the channel rate in the high-SNR regime [7]; such parameters 
are the high-SNR slope (also referred to as the "multiplexing gain") 



<Soo = lim 



Cap(p) 



p— >oo 



\og(Kp) 



and the high-SNR power offset 



lim 



1 



i„ g (*„) - 

P^oolog2 V Soc J 



yielding the following affine capacity approximation 

Cap(p) « ^°° l0g2 (Kp\ dB - 3| d Es£oo) . 

"j|dB 

A direct consequence of Theorem [2] is the following high-SNR characterization. 



Corollary 3 Assume (U$, (H^> and (Wfy. Then = 1 and 

1 



log 2 



1orK-E„ lot 



CoC]|-^7(iV(A = 0)) 
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IV. Reformulations 

We now derive alternative formulations for 'y(N) in Subsection HV-Al and for 7(iV(A = 0)) 
(which characterizes the hign-SNR regime), in Subsection HV-Bl 



A. Non- asymptotic results 

In order to study 'y(N), we express it as the Lyapunov exponent of simpler matrices. For 
i > d + 1, we define the following random matrices. 

( 10 



\ 





Ci , i — d 





1 



Ci,i+d—l 



where Q t i is the coefficient in position (i, I) in Gdn, and set 

d 

Xli = /\ttli. 

Note that (xni)i>d+i an d {fk)i>d+i depend on A. We get the following proposition, whose 
proof is given in Appendix IB. II 

Proposition 4 Assume (HQ), (HE) and Then, N { = n id ■ ■ ■ n^-i)d+i- Therefore, for 

every A > 0, j(N) = d r y(n), hence, 



Cap(p) = log p + E n0t7Td log 



CoC] 



7(n) 



Note that for a given i G N, Ni depends on • • • , C (d+1) \ that is, the fading 

coefficients of 2c? different cells. We now want to reduce the product of the Ni to a product 
of random matrices (that we denote by Hj) depending on the fading coefficients of only d 
cells. Then we reduce it further to a product of random matrices (that we denote by £j) 
depending on the fading coefficients of only one cell. 

By doing so, we achieve two goals: first, we express 'y(N) as the Lyapunov exponent of 
simpler matrices. Second, if the fading coefficients are i.i.d for different cells, then the Sj 
and the are i.i.d. Products of i.i.d random matrices have been studied extensively (see for 
example [8]), moreover, their study can be reduced to the study of a Markov chain on an 
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appropriate space, which can lead to actual analytic expressions (see [5] for an example of 
study of such a Markov chain). 

For i G N*, we denote by Aj the following matrix 







(CiA)- lf 


-QDi + 








{QDi)- 


it 


-{dDiY^DlDi 



and define Hj = f\ d Aj. 

For % > d + 1 , we denote by <5j the following matrix 

( Ci — 



Ci-l,i 
Ci — d,i 



A+|Ci,i| 2 



i.i 



Id d -i 



Oi 



d-i 



d _i 



d-1 



0i„ 



d - 



l,d 



<r~ 

\ ^>i — d,i \ 1, 



d - 



i.i 



^>i — d,i 



CoCi 



and define & = Note that (Aj) i6N *, (Ej) i6N », (<5j)j G N* and (&)ieN* depend on A. 

Proposition 5 Assume (H^j, (HB\) and (Ifify. 

1. For even/ A > 0, 7(H) = j(N), hence, 

Cap(p) = logp + log 

2. ^ssvme K = 1. Then, 

Ai = Sid - • • <J(t_i)rf+i. 
Therefore, for every A > 0, 7 (AT) = d7(£) } hence, 

Cap(p) = log p + E TO log CoCi + 7 (0 
Remark 6 iVote tfia£ for K = 1, for all i eW, 

caAr 1 * = c-^d-v, 
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A; 



Proof: [Proof of Proposition [5] Let us start by proving point 1. We define for i eW 



■a a -acj-Aidd 



Id, 



and 



d ld d 

-it /An rif nt. 



jaA)" T - (CiDi)~ 1 DjDi 

so that for all i e N*, M = P 2 (z + l)A(i)- For * e W ' A i is defined so that A, = P 1 (i)P 2 (i). 
Then, for all n G N*, 

M n . . . Mi = P 2 (n + l)P 1 (n)P 2 (n)Pi(n - 1) • • • P 2 (2)Pi(l) 
= P 2 (n + l)A ri ---A 2 P 1 (l). 



and 



||iV n ...iVi| 



/yp^n+i^-.-^/yp^i) 



< 



A p 2(^+i) 



And) 



Therefore, 7(A r ) < 7(H). Since Pi(l) and P 2 (n + 1) are invertible, we get the opposite 
inequality and point 1 is proved. 

The proof of point 2 is postponed to Appendix IB. 21 ■ 



B. Results in high-SNR regime 

Proposition 7 Assume (HJty, (BJ$) and (Hjjj). Assume moreover that K = 1. 
I. Forie W, we set #J = QD^ and ^ = C~ lf A. Then 



Cr 



log 2 



' ' d—i 



E T log 



CoC] 



- max 

a o<i<d 



7 A* 1 +7 A* 5 
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2. For % > d + 1, we set 



Si— 








/ 


Ci-M 




\ 


Ci — d,i 


Idd_i 




and t/? 2 = 




Ci,i 






Ci,i 










Ci — d , 2 






^ Ci — d,i 


0i,d-i 






1 


Si, it 


M -i 


/ 



Then, 



log 2 



' d—i 



E T log 



Remark 8 1. Recall that for a stationary ergodic sequence of complex random matrices 
(Xi) im * of size d, 



7 



f\x) = and 7 ( /\x \ = ElogldetXi 



2. Note that if for < % < d, TTi = TTd-i and the vectors (O)ieN* are i.i.d, then {ip}) 



i>d+l 



and ((V' J 2 )^)i>d+i have the same distribution and (Vf)i>d+i ((V\ 2 )^)i><m ^o^e i/ie 



same Lyapunov exponents. Therefore, as p goes to infinity, 



log 2 



E T log 



CoC] 



\d/2\ 



ld/2\ 



7 A^' +7 A^ 1 



Proof: [Proof of Proposition [7] Using Corollary [3] and Proposition [U in order to prove 
point 1, we only have to prove that 



7(S (A=0)) = m ? ^ 7 (A*' 
Recall that 7(H) — 71(A) H h 7d(A) and that 

Aj(A = 0) = 



7 



(9) 







(CiA)" 1 * 





By Proposition [2U the sequence 7i(A(A = 0)), . . . , 7 2 d(A(A = 0)) is equal up to the order 
to the sequence 

7 1 (* 1 ),...,7 (2 (* 1 ),7 1 (* 2 ),...,7 d (* 2 ). 
Therefore, © is a direct consequence of 7(S(A = 0)) = 7i(A(A = 0)) H h 7d(A(A = 0)) 



and 7 (A** 1,2 ) =7i(* 



fl,2> 



7 4 (* 
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The proof of point 2 goes along the same lines using the fact that 

Si(\ = 0) 

where 





0d,d 




-ft 







1 



••' 

1 0J 







A 



■• 

1 o oy 



and ^ 





0<2-l,(i-l 




0i, d -i 





therefore, the Lyapunov exponents of ^ 2 and ip 2 are the same. ■ 

V. Bounds on the capacity 

A. Bounds on the top Lyapunov exponent 

We use the Frobenius norm on the matrices, it is a sub-multiplicative norm, therefore, we 
can apply (122]) to the different formulations of the capacity to get the following proposition. 



Proposition 10 Assume (LUty, (hW\) and (fflfy. 
1 . For A = 1 / p and p > 1 , 



Cap(p) < logp + E n(hnd log |CoCj| + ^Elog || jV p (A) • • -iVi(A)|| 



2. For X = 1/p and p > 1 , 

Cap(p) < logp + E 7r0)7rd log 

3. For X = 1/p and p > 1 , 

Cap(p) < logp + E^ 0)7rd log 

4. Assume K = 1 . For X = 1/p and p > 1 , 



CoC 



CoC] 



+ -Elog||tip(A)---ni(A)||. 

P 



+ — Elog||S p (A)---S 1 (A)| 



Cap(p) < logp + log |CoCj| + ^Elog ||£ p (A) • • -6(A)|| . 

Moreover, the bounds are tight as p goes to infinity. 

The bound of point 2 with p = 1 can be reformulated as follows. 

Corollary 11 For X = 1/p, 
Cap(p) < logp+ 
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-Elog 



d- 



2<l 



2d — 2 \ I x 



{l-l)t\d 

Cd+i,d+t+iCi,d +t +i + Ai(i =(J+ i) 

t=(Z-d-l)V0 



+ 



( d ) ( ^ d+1 ' ci+1 ^M+i + Co!+i,2d+iC:L+i,2d+i J 
The proof is postponed to Appendix IB. 31 



B. Other bounds 

Proposition 12 Assume (B\^ and (B\^). For A = l/p, 

max (E, log (A + |Co| 2 ) , K d log (A + |G| 2 )) < Cap(p) - logp < 

Elog(A + |Ci,i| 2 + --- + |C M+ i| 2 ). 

Proof: The upper bound is a consequence of Hadamard's inequality for semi-positive 
definite Hermitian matrices. Indeed, 

1 1 m 

-logdetG m (A) < - ]Tlog(A + |G/ + --- + |CM+/)- 

Ti 777, < * v 



III 



i=l 



Let us show the lower bound of point [12] using the tools of [9] . 



Cap m (p) = — I {x,y\(Ci,i)l<i< m , ■ ■ ■ , ((i,i+d)l<i<m) 

m 



= — y~]l(xj,y \ (arQ !<«,-, (C*,i)i<i< 

m 

> — ^2 1 ( X j'yj-d\( x i)l<i<ji (Ct.i) 
771 . 

m 



mi • • • i 



(Ci,i+d)l<i<m) 
KKmi • • • j (Ci,i+<i) l<i<m) 



which is the per-cell sum-rate capacity of a single user fading channel. Therefore, the lower 
bound is E^ log (l + p |0| 2 ) . 

The role of the distributions ir and tt^ can be exchanged by a right-left reflection, namely 
the transformation Qj = C^m-i+^m+d-j+i- Thus, we get the lower bound. ■ 

In the end of this section, we slightly modify the setting, by considering m cells with 
K single antenna users per cell and m + d single antenna BSts. The communication is 
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characterized by the following (m + d) x Km channel transfer matrix H m , which is a d+ 1 
block diagonal matrix defined by, 

/ Ci,i o ••• o \ 

S2,l C2,2 : 

: : '•• 

Cd+1,1 Cd+1,2 Cm,m 

Cd+2,2 Cm+l,m 



H„ 



\ • • ■ Cd+m,m / 

where are 1 x A' row vectors. 

We consider the per-cell sum-rate capacity Cap m (p) that is given by (CQ). We denote by 
Cap(p) the limit of Cap m (p) as m tends to infinity. 

Note that in the limit, this setting is equivalent to the setting we define in Section [Til In 
particular, the normalization by 1/m or l/(m + d) is equivalent. 

Proposition 13 For all n £ N* ; 

Cap n ( — — p < Cap(p) < Cap n (p). 



n + d \ n 
Moreover, the bounds are tight as n goes to infinity 

Note that taking the upper bound for n = 1, one gets the upper bound of Proposition! 
Proof: Number the cells from 1 to m and the antennas from 1 to m + d. 

Upper bound. Take 1 < n < m. Consider the following transformation of the communica- 
tion system: duplicate the antennas number n + 1 to n + d (each antenna is replaced by two 
antennas) and denote by (n + 1, a), (n + 1, b), . . . , (n + d, a), (n + d, b) the new antennas. The 
cells number n—d+1 to n broadcast toward the antennas (n+1, a), (n+2, a), . . . , (n+d, a) and 
the cells number n+1 to n+d broadcast toward the antennas (n+1, b), (n+2, b), . . . , (n+d, b). 
See Figured] for an illustration with d — 2, m = 8 and n = 4. The new system has a higher 
capacity since a first step of the decoding could be summing up the signals received at 
antennas (n + i, a) and (n + i,b) (1 < i < d) to get the signal received at antenna n + i in 
the former system. Therefore, 

mCap m (p) < nCap n (p) + (m - n)Cap m - n (p). 
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Take k G N*, by induction, 

(nk)Cap nk (p) < (nk)Cap n (p). 

Dividing by nk and taking k to infinity in the LHS gives the upper bound. 

Lower bound. Take k G N* and consider k(n + d) users and their corresponding antennas. 
For every group of n + d users, silence the d last users and redistribute their power to the 
n first users so that their SNR becomes (n + d/n)p (the average SNR is still p). See Figure 
[2] for an illustration. Since asymptotically, the equal power distribution among the users is 
optimal [10, Appendix C], we get that the new system has a lower capacity. Therefore, for 
k going to infinity, 

(n + d \ 
— — p J < k(n + d)Cap k{n+d) (p). 

Dividing by k(n + d) and taking k to infinity in the LHS gives the lower bound. ■ 

C. Numerical comparison of the bounds 

We first compare the bound of Corollary [11] and the bounds of Proposition [121 Note that 
by the ergodic theorem, the upper-bound of Proposition [12] grows like log d, whereas in 

2 

the bound of Corollary [11], the part log ( 2d ^ 1 ) Cf+i,d+iCi d+i alone already grows like d. 
Nevertheless, it is not necessarily true that the upper-bound of Proposition [12] is better that 
the one of Corollary [TT] for all d and all fading distributions. 

In Figure [3] we present the bounds of Corollary [11] and Proposition [12] in the special case 
of Rayleigh fading (real and imaginary parts are independent Gaussian random variables 
with zero mean and variance l/y/2). The curves are produced by Monte Carlo simulation 
with 10 6 samples. We see that in this case, even for d small, the upper-bound of Proposition 
is better than the one of Corollary [11] 

In Figures [4] and [5] we compare the bounds of Proposition [TO] point [T] and [3] Proposition 
and Proposition [13] in the special case of Rayleigh fading (real and imaginary parts are 
independent Gaussian random variables with zero mean and variance l/y/2). The curves are 
produced by Monte Carlo simulations with 10 5 samples. 

Note that in the case d = 2, for K — 1, the bounds of Proposition [TO] point [3] are better 
than those of point [1] whereas for K > 1, it is the opposite. 
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We see that in the case d = 2, for K = 4 and K = 10, the upper-bound of Proposition 
[121 is very close to the capacity and the upper-bounds of Proposition [TOim are getting tight 
very rapidly. 

In the case d = 2, we want to compare the random-fading channel with the non-fading 
channel. See Appendix [D] for the per-cell sum-rate capacity of the non-fading channel. The 
comparison is done in Figure [5} in the eight cases that we consider, the random-fading 
channel is better than the non random one. 

VI. Results for particular cases in the high-SNR regime 

A. Case d = 1 

As a direct application of Proposition [7] in the case d = 1 and K — 1, we get the following 
result. 

Proposition 14 Assume (!{J\), (B~E) and (B~M). Then 

too = t 2 max (E no log | Co I 5 E*i log |Ci|)] . 

Note that a similar result was already proved by other techniques in [5] under much 
stronger hypothesis, in particular, independence of the fading coefficients was assumed there. 
In contrary, our result depends only on the marginal distributions of the fading coefficients 
and is valid for a larger class of joint distributions. 

We want to compare the per-cell sum-rate capacity of the random-fading and non-fading 
channels. For a random variable £, by Jensen's inequality, 

ElogKI 2 < logEKI 2 . 

Therefore, under the constraints E„. |£o| 2 < 1 and |d| 2 < 1, the non-fading channel 
achieves the best per-cell sum-rate capacity in the high SNR regime. 

B. Case d = 2 

We now assume that d = 2 and K — 1 and that the fading coefficients have the following 
form; for % e N*, 

Ci-2,i = aat , Ci— i,i = and Q,i = Ci, 
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where aj, 6j and Cj are random variable distributed according to 7r a , 7r^ and tt c respectively 
and a and (3 are parameters such that a > and /? > 0. Moreover, take the following 
normalization that can always be achieved by modifying a and (3. 

E^loglail = E^logl&il = E^loglcil . 

We use the notation of Proposition [7J 

Proposition 15 Assume that (a^, bi, c.j)j e N* a stationary ergodic sequence such that for 
all i G N*, almost surely, ai and Cj are non zero and #ia£ t/iezr exzsi e > such that 
E^ (log K|) 1+£ , E^ (log |foi|) 1+e and E Wc (log |ci|) 1+£ are finite. 

Then, there exist a domain V C (0, 1] x [0, 1] such that for all (x, y) G V, (0, x) x [0, y) C V 
and for all (a, (3) G T>, as p goes to infinity, 

—2 

£co = ] — -E* a log|a|. (16) 
log 2 

The proof is postponed to Appendix IB.41 

Remark 17 1. The set T> is not maximal in the sense that (To)) may hold for couples 
{a,P)£V. 

2. Note that for (a, (3) G T>, in the high-SNR regime, the lower bound of Proposition [7H 
is tight. 

3. The proof will yield an effective construction ofD, which allows us to find many points 
in T>. Indeed, we construct (/ p ) p gN* o, family of functions on (0,1] x [0,1] with the 
following property: if there exists p6N* such that f p (a,/3) < 0, then (To)) holds. 

4. Note that [To]) does not hold when a > 1, indeed, as it will appear in the course of the 
proof, as p goes to infinity, 

Cap(p) > logp + 2E 7ra log \a\ + log a. 

We conjecture that [To]) does not hold when (3 > 1 either. 

Let us apply Proposition [TBI to the case where (aj, foj, Ci)^* are independent Rayleigh 
distributed coefficients. In Figure El we plot points for which / 2 o is less or equal to -0.05 
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(Monte Carlo simulations realized with 10 5 samples). Therefore, ffl6|) holds for (a, (3) in the 
stripped region and in particular for a, ft < 0.4. Note that in this case, the power offset is 

7 



r _ 

log 2 



where 7 is the Euler constant. 



C. Artificial fading 

In the frame of non-fading channels, we consider artificial fading, that is, every user uses a 
pseudo-random fading and multiplies its signal by this artificial fading. The fading coefficients 
then have the following form, for i e N* and < s < d, 

where cko, • • • are non random positive numbers and Pi , i G N* are stationary ergodic 
pseudo-random complex row vectors of size K distributed according to a law denoted by 
7Tp. We moreover assume that for all i £ W, almost surely, the coefficients of Pi are non 

1 1 2 

zero and that E„. p ||Pi|| =1. 

In [11], it is proved that in the case d = 2, the per-cell sum-rate capacity is smaller with 
artificial fading. Indeed, had such a procedure helped, then it would be used in non-fading 
situations to enhance capacity. It is evident then that it is deleterious, as the expression in 
Proposition [18] exhibits. 

We consider the high-SNR regime and derive the explicit influence of the artificial fading. 

Proposition 18 Denote by the power off-set without artificial fading (that is, Pi = 
(1, . . . , 1) almost surely) and by the power off-set with artificial fading. Then, 

C oo = C oo -^—E 7Tp \og\\P 1 \\ 2 . 
log 2 

Remark 19 By Jensen's inequality, we get that > C ^, therefore, in the high-SNR 
regime, the per-cell sum-rate capacity is smaller with artificial fading. 

Proof: [Proof of Proposition [18] We set until the end of the proof A = 0. Using Corollary 
[3] and Proposition [SJ 



C 



00 



log 2 



logK -1L lot 



a a\ 



log II-P1H 2 - 
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whereas 



log 2 



\ozK-E„ lot 



a a d 



denote the matrices without artificial fading. Therefore, we only have to 



where 

prove that for z G N*, A, does not depend on (Pj)jgN*. In the case if = 1, Aj = 5 id ■ ■ ■ %-i)d+i 
and for z > d+ 1, 5* does not depend on (Pj)j £ N*, therefore, Aj does not depend on (Pj)j 6 N*- 
Let us assume K > 1. Using Proposition [5l we only have to prove that for z G N*, Aj 
does not depend on (Pj)^*. 



and 



aoPd(i-l)+l a lPd(i-l)+2 • • • Oi d _iP d ^ 
a Pd(j-l)+2 ' ' ' &d-2Pdi 

••• a Pdi j 



^ ®dPd{i-l)+^ ^d-lPdii-^+l^ ■ ■ ■ OCiPd(i-l)+l 



Di 



«dPd(i-i)- 



«2Pd(i-l)- 



V 



a^Pdi 1 " y 



Let us define another channel transfer matrix H m by K = 1 and for z G N* and < s < d, 

Ci,i+s Q^s l|Pi+s|| ■ 

In the same manner, we define Cj, -Dj and Aj. A straight forward verification shows that for 
z G N* 

<7iC} = CtDi = CiDi and D+A = 5jA. 

Moreover, since Aj is a function of CjCj, Cj-Dj and D\D^ Aj = Aj. However, since if = 1, 
we have already proved that Aj does not depend on {Pi)%eWi therefore, Aj does not depend 
on (Pj)j eN ». ■ 
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VII. Numerical simulations 

A. Influence of the correlation 

We assume that the fading coefficients are Rayleigh distributed (real and imaginary 
parts are independent Gaussian random variables with zero mean and variance l/y/2) and 
independent for different users. We are interested in the following question, which of the non- 
fading channel and the Rayleigh fading channel gives a higher per-cell sum-rate capacity. 

In the case d = 2, A = 0.1, 1 and K — 1, 2, 4, 10, with all fading coefficients independent, 
it is shown in Subsection IV- CI that the Rayleigh fading is beneficial. 

In the case d = 1 if we assume independence between the Qj, it is known that Rayleigh 
fading is beneficial over non- fading channels in the high-SNR region already for K — 2 ([5]). 
If we assume that for % G N*, d-i,i = then, the sum-rate per-cell capacity is less than 
the one of a non-fading channel (see Subsection IVI-CI and [11]). We investigate the following 
question: what is the maximal level of correlation between £i-i,i an d Cm that still provides 
benefit over the non-fading channel. See Appendix [D] for the derivation of the capacity of 
the non-fading channels. We denote by c the correlation between the real (resp. imaginary) 
part of and the real (resp. imaginary) part of Ci-i,i 

In Figure [7] we present the bounds of Proposition [iulfTl and Proposition [13] in the special 
case of Rayleigh fading. In Figure M we present the bounds of Proposition [TOTTI and Propo- 
sition [13] in the following special case: is Rayleigh distributed, is a G [0, 1] times 
a Rayleigh distributed random variable. In both cases, the curves are produced by Monte 
Carlo simulation with 10 5 samples. 

We see that even with a correlation close to 1, fading still provides an advantage over 
non-fading channel. Moreover, note that K large, high SNR and a close to 1 are conditions 
in which the advantage of the fading is larger. 

B. The asymmetric Wyner model 

With the following specification, the model studied is the Rayleigh-fading Wyner model 
([3]). We take d = 2 and the Qj independent with the following distributions. For i G W, 
d t i+i is Rayleigh distributed (real and imaginary parts are independent Gaussian random 
variables with zero mean and variance and £t,i (resp. Q,i+2) is a G [0,1] times a 
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Rayleigh distributed random variable. The asymmetric (Rayleigh-fading) Wyner model is 
similar to Rayleigh-fading Wyner with a slight modification. For i 6 N*, is Rayleigh 
distributed and Cm+i (resp. Q,i+2) is a times a Rayleigh distributed random variable. Note 
that in Subsection IVI-BI we prove that in the asymmetric case, the power offset for a < 0.4 
is 7/ log 2. 

The two models are very similar and yet, in the non-fading case, the per-cell sum-rate 
capacity is notably different (see Appendix [D] for the derivation of the capacity of the non- 
fading channels). In Figure [9] we present the capacity of the two models without fading and 
the bounds of Proposition [13] for the two models with Rayleigh fading. We study one case 
in moderate SNR (A = 1) and one case in high SNR (A = 10~ 4 ). The curves are produced 
by Monte Carlo simulation with 10 5 samples. 

Note that in the high-SNR region, for the non-fading channel, the per-cell sum-rate 
capacity is very different for symmetric and the asymmetric models, whereas the per-cell 
sum-rate capacities for the symmetric and asymmetric Rayleigh-fading models are very close 
(but not equal as shown in Figure [TOl for A = 10~ 4 and a = 0.5). 

To understand better the influence of fading on the difference between the two models, 
we present in Figure [11] the bounds of Proposition [13] for the capacity of the two models 
(symmetric and asymmetric) with the following fading: the modulus is uniformly distributed 
between 1 — e and 1 + 5 and the phase is uniformly distributed between and 2ett, where e is 
a parameter between and 1. Note that for e = 0, there is no fading and for e = 1, the fading 
is uniformly distributed on the disc of center and of radius 2. The curves are produced by 
Monte Carlo simulation with 10 5 samples. We notice that the difference between the two 
models decreases between e = and e = 0.5 and that in high-SNR, it increases slightly 
between e = 0.5 and e — 1. 

VIII. Concluding Remarks 

In this paper, we study the per-cell sum-rate capacity of a channel communication with 
multiple cell processing. The main tools is a version of the Thouless formula for the strip 
which we prove in the article. It allows us to prove that the per-cell sum-rate capacity 
converges as the number of cells and antennas goes to infinity. We give several expressions 
of the limiting capacity in terms of Lyapunov exponents and several bounds on the per-cell 
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sum-rate capacity. 

We apply those results to several examples of communication channels and get insight 
on the evolution of the capacity as a function of the key parameters of the problem. In 
particular, in the high-SNR regime, some explicit formulas are derived. 

Note that the model here applies verbatim to randomly varying intersymbol interference 
channels. 

Some of the tools of this article can be used to derive CLT-type results on the capacity 
in order to study the outage-probability. Details will appear elsewhere [12]. 
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A. Random Schrddinger operators techniques 

1 ) Lyapunov exponents theory: We use the theory of product of random matrices. For a 
general introduction to the aspects of the theory we use here, the reader may consult [8], 
[13], [14], [15], [16] or [17]. See appendixIClfor the relevant background on exterior products. 

Theorem 20 (Furstenberg H., Kesten H. (I960)) Consider a stationary ergodic se- 
quence of complex random matrices (Xi)^* of size p and any norm on the matrices. Assume 
moreover that 



Appendix 



Elog + ||Xi|| < oo 



then a.s, n 1 log ||X„ ■ • • Xi\\ converges to a constant: 



lim -]og\\X n ..-X4±>y{X). 



We define p constants 71 (X), . . . , 7 P (X) such that for 1 < i < p, 
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Proposition 21 

7iP0>--->7 P P0- 

The constants 71 (X) > ■ ■ - > 7 P (X) are called the Lyapunov exponents and ^(X) = ^yi(X) 
is called the top Lyapunov exponent. 

We will also use the three following properties: 

1. For any sub- multiplicative norm, for p6 N* 

7 (X) < ^E\og\\X p ---X 1 \\, (22) 
and the limit of the RHS as p goes to infinity is 'j(X). 

2. 

-ElogldetXJ < 7 (X). (23) 
p 

3. Assume that the matrices (Xi)^* are i.i.d, then for all 1 < i < p, ji(X) = ji(X<). 
Finally, we quote the following proposition [18, Proposition 1]. 

Proposition 24 Consider a stationary ergodic sequence of complex random matrices (Xi)^* 
of size p and any norm on the matrices. Assume moreover that 

Elog + IIXJI < 00. 



2^ 

igN* ; 



Finally, assume that there exist three sequences of random matrices (X^)^*, (X, 
(Xf)i(zF}*, of respective sizes k x k, (p — k) x k and (p — k) x (p — k), for 1 < k < p — 1, such 
that almost surely, for all i 



X, 



x} 






xf 



Then, 71 {X), . . . , 7 p (X) is equal up to the order to the sequence 

7l (X 1 ),..., 7fc (X 1 ),7i(X 3 ),..., 7p _ fc (X 3 ). 
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2) Proof of Theorem^ 1 : In order to prove point 1 of Theorem [2], we first prove a slightly 
more general lemma. 



Lemma 25 Assume (1%]$), (ffify and (B^j. For all X e C such that X £ 



almost surely, 



— log|det G dn \ 



- 3Elog|det(C 2J D 2 )| + 37 (N) . 



dn n-*co d d 

Proof: For i e W, set B, t = QCj + £>J +1 A+i + Xld d and A { = C.D,,. Note that the 
eigenvalues of Gd n are bounded away from zero. To compute log |det Gd n \, we write the 
following decomposition: GdJJdn = L dn , where Ud n is the upper triangular by block matrix 



/ 



dn 



x 1 


x 1 




x 1 


o d 


x 2 




x 2 










Od 




o d 


x n 



\ 



the Xi are d x d matrices such that X = d , X\ = Idd, and for % > 1, 

MXi-\ + BiXi + A\ +l X i+ i = d . 
L dn is the lower triangular by block matrix 

\ 



(26) 



-a\x 2 


o d 




o d 


A 2 X 1 


-4x 3 






o d 






o d 


o d 


o d 


A n X n _i 


~ A n +lX n +l 



That decomposition allows us to write log |det Gd n | as a determinant by block, 

n n n 

log | det G dn | = y^log) det A i+1 | + ^ log | det X i+l \ - ^ log | det X { \ 



i=i 

n 



^ log | det + log | det X n+X \ . 



i=i 



Therefore 



1 1 n 1 

— log | det G dn \ = — ^ log | det A i+1 1 + — log |det X n+1 \ . 
dn dn dn 



(27) 



X^iLil°g | det Ai+i | converges by ergodicity toward E log |det A 2 \. Note that the choice of 
A n+ i is arbitrary, indeed, if we take another value, say A n+ i, then A\ +1 X n+ i = A\ +l X n+ i 
and ( |27|) stays unchanged. 
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We emphasize that the derivation of (|27j) is inspired by Narula's thesis ([19]). 
The Xi are defined by (1261) . We can reformulate it in the following way. Set V$ = 
then (I26I) is equivalent to V$+i = M^/i and moreover, 

d 



Xi-i 
X, 



X n 



+1 



0, 



IcL )M n ---M 1 



Id, 



Denote f — f\ 
O We get 



0d 



. For the relevant background on exterior products, see Appendix 



AyV n _ 1 ...iV 1 (i; 1 1 A---A<). 
However, v\ A • • • A vf = f. Therefore, 



a a I 

detX n+1 = /\X n+1 = J\ M Q c 



Id d )M n ---M 1 



Id 



(2« 



= / t iV n iV n _ 1 ...iV 1 /. 

Taking the canonical basis of f\ d C 2d , / is the last vector of the basis and det X n+1 grows 
like the bottom-right coefficient of the product of the N iy therefore, its growth rate is bounded 
above by the Lyapunov exponent of the iVj. 



limsup-Elog|detX n+1 | < 7(JV). 



(29) 



Using (1271) . it is enough to prove the opposite inequality to conclude the proof. The end of 
the proof is inspired by [20]. 

Lemma 30 If there exist a basis of /\ d C 2d , say (gi)iei, such that for all i,j G /, almost 
surely, 



lim inf — log 

n— >oo 77, 



g]N n ---N igi 



< lim inf — log [det X n+ i \ 

rwoo 77, 



(31) 



then, almost surely, 



7 (A^) < lim inf - log [det X n+1 \ . 

n—*oo 77, 



Let us first prove the lemma. 
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Proof: For any finite basis Si and £2 in a vector space, we have for all A 

sup WA/3\ > c\\A\\ 
aeS!,pes 2 

for some universal c. Thus, ( 13T1) shows that, almost surely, 

j(N) < lim inf - log |det X n+1 \ . 



To finish the proof of Lemma [25], we denote by {ei, . . . , e 2d } the canonical basis of C 2d 
and we apply the lemma with the following spanning system of f\ d C 2d 

S = {(ei + ei) A ■ • • A (e d + e d ) , e 1 , . . . ,e d e vect(e d+ i, ■ ■ ■ ,e 2d )}. 

For a choice of ef , . . . , ej, such that for 1 < j < d, ej = 2Ji=i a ij' e d+i) we define the 
d x d matrix of the a^j. We get 

gi = (e 1 + et)A---A(e d + e*) = /\ 



Ma 



In the same way, for a choice of e\ . . . , e b d e vect(e<i + i, . . . , e2d), we define a dx d matrix, 
such that 



g 2 ±(e 1 + e b 1 )A---A(e d + e b d )=f\ 
We define two new sequences (AA and (J5 4 ) such that 



E h 



. For 2 < i < n- 1, B i = B i; 
• For 1 < z < n, Aj = Aj, 



. Bi = -A 2 £# 
• -Bn = A n (E b ) 
. A 



igN* 



fN n N n _ 1 ...N 1 f = ff\ 



o d 


ld d 


ld d 





iGN* 

</ 



r ]N n ^...N 2 /\ 



Then, 




o d 


Id, 


-TiT" 11 " 


E* 



f 



giN n _ x . . . N 29l . 
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Therefore, to prove the condition (13T1) . it is enough to prove that, almost surely, 



lim sup — log 

n 



detG 



(in 



-log|det G dn 

n 



< 0. 



We now use perturbation theory techniques. Indeed, we denote by p the spectral radius of 
a matrix, i.e. its largest eigenvalue in absolute value. Recall that for a matrix S, p(S) < 
y/p(SS*) and that y/p(SS*) is a sub-multiplicative norm. As a consequence, for positive 
Hermitian matrices, the spectral radius is sub-multiplicative. Moreover, we denote by 



the Frobenius norm. Recall that y/p(SS^) < \\S\\ F . We will also use the fact that the 
eigenvalues of G dn are bounded away from 0by/i = AifA>0oryU = if A ^ K. 
Moreover, we define U dn = Gd n ~ Gdn-, which has rank less than or equal to 2d. 



-log 

n 



det G dn 



-log|det G dn \ 
n 



- log |det(G dn + U dn ) \ - - log |det G dn \ 
n n 

-logldet^d^+G- 1 ^)! . 
n 1 1 



G d ^Udn has rank at most 2d, therefore, 



n 



log 



detG 



dn 



1 2d 

- log | det G dn \ < — log 1 + p(G^U dn ) 
n ti 



2c/, 

< — log 

n 

2d 

< — log 

n 

2d, 

< — log 

n 

2d, 

< — log 

n 



1 
1 
1 

1 + 



P( G dn U dnU dn G Q 11 



dn 



P(G- d ^G- d ^p{U d nU t 



t ^ 

dm 



-.^PiUdnUl) 



\u, 



dn || p 



Moreover, 



\TT II 2 

I u dn 1 1 p 



Si + Xld d +t\e* 



hence, with the integrability condition, sup n Elog 
inequality, for a given r\ > 0, 



P 



1 



\S n + Xld d +TlE 

l+e 



b\\ 2 
F 



U, 



dn\\ p 



< oo. By Tchebicheff 









( - log 


1 H \\Udn\\ p 


> v) 




P 



sup n Elog 


1 + ji \\U dn \\p 


l+e 







The RHS is a summable series, therefore, by Borel-Cantelli Lemma, almost surely, 



lim sup ( — log 

n 



detG 



dn 



-log | det G dn 

n 



< 0. 
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This finishes the proof of Lemma [25j ■ 

Remark 32 By the same kind of perturbation theory techniques, we can show that in order 
to study the limit in m of Cap m (p), it is enough to study the sequence every d steps. 

For a hermitian matrix h whose ordered eigenvalues are a%, . . . ,a n , we denote by the 
spectral distribution of h, the measure 




i=l 

where 5 X is a Dirac measure at x. 

The following technical lemma will be used several times to prove domination properties. 



Lemma 33 Denote by \i n the spectral distribution of H^ n W dn . Consider the following diag- 
onal by blocks matrix: 



dn 





( 2£i 


o d 




o d 






2B 2 














o d 




\ o d 




o d 


2B n 



and denote by fi n its spectral distribution. Then, for any non-decreasing function f, 

fd\i n < \ fdfi n . 



Proof: Denote 



H, 



tin 





-D\ 


O d ,dK 




Od,dK 


Od,dK 


C 2 


~D\ 














Od,dK 


d}d K 




O d ,dK 


C n 





\ 



then F dn = H dn H dn + H dn H dn . Since Hd n H\ n is a non-negative Hermitian matrix, by Weyl's 
inequalities, for all 1 < i < dn, the i-th eigenvalue of Hd n H dn is less or equal than the i-th 
eigenvalue of F dn - ■ 
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CoCj 



From Lemma 



First note that (l/d)Elog |det(C 2 D 2 )| = E^log 
A > 0, 

1 1 n 

— log | det G dn \ < log 2 + - V log | det B { 



we deduce that for 



n 



i=l 



Therefore by (H2J) arid Hadamard's inequality, (1/dn) log|detGdn| is a uniformly integrable 
sequence and the almost sure convergence of Lemma [25] implies point 1 of Theorem [2l 
3) Proof of Theorem^ 2: We begin by a few notations. For A ^ R, set 



/(A) 



lim — lo£ 

n^oo an 



det 



Aid 



c/fi 



which exists by Lemma [25l The existence of the weak limit of \x n and the fact that it is non 
random is a classical fact of the random Schrodinger operators theory, see for example [16, 
Theorem 4.4]. For A G C, we set (if it exists) 



<?(A) 



log | a; + A| d\i{x). 



We emphasize that since log is not a bounded function, we cannot directly deduce from 
Lemma [251 and the weak convergence of the \i n to fi that for A ^ WL, /(A) = g(\). 
Finally, for A G C, define 

/i(A)=E.log|CoC]| + ^7(AO- 
The following lemma is a generalization of the Thouless formula for the strip proved in [20]. 



Lemma 34 Assume (I^j, (Iflfl and For all A G C, g(X) = h(\). 

The proof of this result is done in the frame of channel transfer matrices but one does not 
need to assume that the Ai and the Bi are upper triangular by blocks, one just need instead 
of (H[3]) the hypothesis that almost surely, AiBi is invertible. 

Proof: The proof goes along the following lines, we first prove that for A ^ R, g(X) 
exists and equals to h(X), then, following [21] we argue that g and h are two subharmonic 
functions equal everywhere except a set of measure, therefore they are equal everywhere. 

Step 1: Let us first prove that for A ^ M, g(X) is well defined. log|x + A| is bounded 
away from — oo, therefore, g(X) exists although it may be oo. For R > 0, let us denote 
by log R the function t — > log(t) A R. By monotone convergence, it is enough to prove 
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that Jlog^ \x + A| dfi(x) is bounded uniformly in R. Since x —> log R \x + A| is a bounded 
continuous function, 

/log R \x + A| dfi(x) = lim / \og R \x + X\ dfi n (x). 
By Lemma [331 an d using that J log R \x + A| dfi n (x) < J log |x + A| d[i n (x), 

lim / log^ |x + A| dfi n (x) < lim / log \x + A| djl n (x) = Elog |det£>i| < oo, 

n— >oo / n— *oo / 

where the last inequality comes from (H2D and Hadamard's inequality. Finally, we get that 
for A £ E, 

#(A) < Elog|detBi| < oo. 

S^ep Let us prove that for A ^ R, /(A) = g(A). Applying Lemma [331 one shows that 
for A ^ C, the sequence (/log \x + A| <i/i n (x)) rigN , is uniformly integrable and therefore, 

/(A) = lim E / log \x + A| dfx n (x). 

n— >oo / 

By Lemma [33l for R > 0, 

E / log |x + A| dfi n (x) < E / log \x + A| dfl n (x) 

Jx>R Jx>R 



E 



/ log |x + A I dfii(x). 

Jx>R 



Therefore, for n G N* and R > 0, 



E / log |x + A| dfi n (x) — g(\) < E / log |x + A| d^i{x 

J Jx>R 



+ 



+ 



E J log R \x + A| dfi n (x) - 7, J log R \x + A| dfi(x) 
E / log \x + A| d/ji(x) 

Jx>R 



We first fix R > such that the first and the third terms are arbitrary small and then, by 
weak convergence, the second term goes to as n goes to infinity. Therefore, for A ^ K, 
/(A) = g(X) and by Lemma [251 #(A) = h(\). 

Step 3: Let us prove that g and h are subharmonic on C. See [21] for the relevant 
definitions. Since for i e N*, iVj(A) is an entire function of A, is subharmonic ([21]). 

Let us prove that g is subharmonic. For R > 0, set 



H (A) = y (log|a; + A|V-fl)d/i(ar) 
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By Lemma [331, gn a continuous function. As R goes to infinity, g^ is a decreasing sequence 
of functions converging point wise to g, therefore, g is subharmonic. 

The functions g and h are subharmonic on C and equal on C - R, therefore, g and h are 
equal on C. ■ 

To finish the proof of point 2 of Theorem [21 let us prove that h(X) converges to h(0) 
when A goes to in R+. Note that Elog |det Mi| = 0, therefore, Elog|detiVi| = 0. By 
(E3D, h(0) > E n \og CoC] , therefore, using Lemma [33] and the fact that for X,x £ M + , 
log | a; + A | > log \x\, we get the desired result. 



B. Other proofs 

1) Proof of Proposition^ We use the notation of Subsection IA.2I We define Xj for % £ N 
and 1 < j < d such that the the element at the position (s,t) of Xj is xf 1 - )d+s . Recall 



that GdnUdn — Ldn- Therefore, for a given j such that 1 < < d, we get the following 
characterization of the sequence = for — d+ 1 < i < 0, x* = 5jj for 1 < % < d and 



for i > d + 1 , 



Therefore, 



Moreover, 



E c«4 = o. 

l=i— d 



( x i ~ d + l \ 



TtX; 



( x {i-l)d+l 



(i+l)c! 



V «1 



mid • ■ ■ m(j_i) rf+ i 



(i+l)d 



x 



( (i-2)d+l 



V 



= rriid • • • m.(i_i)d+iVi. 
Therefore, together with [26l it proves Proposition [U 



.r 



(i-2)d+l\ 



.r 



id 



(35) 
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2) Proof of Proposition^ 2: In order to prove point 2 of Proposition [5], we first prove the 
following lemma: 

Lemma 36 For all i > d + 1, there exist matrices p\(i), p|W f or 1 — s — d an d fi s {i) f or 
1 < s < d + 1 such that 5\i) = n(i), 5 s {i) = <f (0, +d ~ s+1 ), p»(i) = . . . , 

and p s 2 {i) — p s (C +d ~ s ), where 5 s , p^ and p s ^ are deterministic functions. We have moreover 
the two relationships 

8 s {i)=p s 2 {i + l)p\{i). (37) 
S s+1 (i)^pl(i)p s 2 (i). (38) 

Finally, for i > d+1, 5i = 5 d+1 (i) 

Proof: For z > d + 1 and 1 < s < d, define 
. for s < I < 2d, 

(l-s)A(d-s) 
t=(l-s-d)V0 

. for 1 < I < d, a\~ s = -tL.+i-i,i+d-./(i+d-sj+d-s> 

• for 1 < / < S, b\ s = — Ci-s+l,i+d-s/ Ci-s,i+d-s, 

• P t/Ci-s,i+d-sCi+d-s,i+d-s- 

Then 





/ 


Id s _i 




s -l,2d-s+l 


\ 


m = 


Oi,,_i 


a 


i s *> s 
s ' ' " a 2d 








2d - 


s,s 


Id 2 d- S 


/ 



and 
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\ 



Oi 



d-i 



a 



Finally, for 2 < s < d, 
/ 



b i-s+l 
yi 8+1 


ld s -2 


s -2,2d-s+l 


02d-s+2,s-2 




i,s 


02d-s,l 


a 2d 


2 d-s,l 


ld 2 d-s 


pi-8+1 


n „ i— s+1 _ i—s+1 

Oi,d- s +i a x ■ ■ ■ a d 



V 

A (straight forward yet tedious) verification shows that ( 1371) and ( l38l) are satisfied. 



Note that in the proof, we make a choice of particular p s , p s 2 and 5 s . Point 2 of Proposition 
|5] is a direct consequence of the following lemma 



Lemma 39 For all i G N*, 



Therefore, 



^■i — b~id - ■ ■ 



— i — iid ■ ■ ■ £(i-l)d+l- 
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Proof: With the matrices of Lemma [36], we can transform the product of the /ij using 
alternatively (157)) and (158)) . 

M i = P 2 (z + l)P 1 (z) 

= P 2 (z + l)A(z) (J*,(i)) -1 , 

= 5 1 (id)---5\(i-l)d+l) 

= p\{id + l)pl(icZ)pi(id)pl(id - 1) ■ • -pi((z - l)d + 2)pJ((z - l)d + 1) 
= p\{id + l)5 2 (zd) ■ ■ • 5 2 {{i - l)d + 1) (p\{{i - l)d + l))" 1 

= p\{id + 1) • • -pg(z<i + l)5 d+1 {td) ■ ■ ■ 5 d+ \{i - l)d + 1) [p\{{i - l)d+ 1) • • -^((i - l)d + 1)] 

where the last equality is proved by induction. 
Therefore 

P 2 (z + l)A(z) (P 2 «)^ = 

p 2 (zrf + 1) • • -^((id + l))6 id ■ ■ ■ « (w)w - l)d + 1) ■ • - l)d + 1))" 1 

and 

[p 2 (*d + 1) • • -p d 2 ((id + 1))] - 1 P 2 (t + 1) = 

(40) 

• • • 5 {i - 1)d+1 [p\{{i - l)d+ 1) • • -^((i - l)d + 1)] P 2 (z) (A(z))" 1 
At this point, we emphasize that their exist a deterministic matrix valued function A such 

that for all i e N*, A t = A (C d(i_1)+ \ • • • , C*)- In the same wa y' we define £l and The 
RHS of (HO)) is a function of . . . ; Q d whereas the LHS is a matrix valued function of 

( id+1 , . . . , thus, both functions are constant. Therefore, there exist a matrix X such 

that for all % 6 N* 

P 2 (* + 1) = p*(id + 1) • • + 1))J. (41) 

Therefore 

A, = T~ x b~ id ■ ■ • 5^i) d+ iX. 
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Note that (jUJ) can be rephrased in the following way. P_ 2 and v\' ' 'V\ are equal up to 
multiplication by a constant to X. Therefore, to prove that X = Id 2 <z for the choice for p s , 
p s 2 and 5 s that we have made in Lemma [36], it is enough to prove that for one given value 

of c\--.,c d , 

p 2 (c\...,c d ) =£(( d )---pX)- 

We will prove it for C 1 = • • • = ( d = (1, 0, . . . , 0, 1). Indeed, 

' d Id, 



P 2 ((1,0,...,0,1),...,(1,0,...,0,1)) 



id, -id. 



For 1< s < d, 





/ 








p*((l,0,...,0,l)) = 


2 <2-l,l 


Id 2 d-i 




( 1 


0i,d-i 


-1 


M -i 



Hence, by induction onKKd, pj((l, 0, . . . , 0, 1)) • • 0, . . . , 0, 1)) 









\ 






Id 2 d-t 




Id* 




-Id t 





V 

Therefore, X = Id 2 d- ■ 

3) Proof of CorollaryXTDi Let us compute E log Hn^+i || . To that extent, we define (ei, . . . , e 2 ,) 
the canonical basis of C 2d and we take (e^ A ■ ■ • A ejjl < i\ < ■ ■ ■ < id < 2d) as a basis 
of C(^). For given 1 < i\ < ■ ■ ■ < i d < 2d and 1 < ji < ■ ■ ■ < j d < 2d the coefficient of 
n<m( e u A • • ■ A e id ) in A • • • A e Jd (we denote by a its absolute value) is the determinant 
of the d x d sub-matrix of Hd+x obtained by taking the lines 1 < ji < ■ ■ ■ < j d < 2d and the 
columns 1 < %\ < ■ ■ ■ < i d < 2d; we denote the latter sub-matrix by T>. Denote by the 
coefficient at position (i, I) of Gd n - 
. If l<jt <■■■ <j d <2d-l, 

— if for all 1 < s < d, i s — j s + 1, then a — 1; 



otherwise, there exists a line of zeros in V, therefore, a = 0. 
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• If 1 < ji < ■ ■ ■ < jd-i <2d-l, and j d = 2d, 

— if there exists 1 < Sq < d — 1 such that for all 1 < s < So, i s = j s + 1, for all 
s < s < d, i s = j s -! + 1 and j SQ = I g {ji, . . .,j d -i}, then a = (d+i,i/(d+i,2d+i ! 

— otherwise, there exists a line of zeros in V, therefore, a = 0. 

We now count how many times each value appears as the absolute value of a coefficient 
of n d+1 . 

• To pick 1, one needs to pick d lines among the first 2d — 1 lines of fi d +i and then, one 
has no choice for the columns: ( 2d d 1 ) choices. 

• To pick Cm,i/Cm,2fM-i , one needs to pick d— 1 lines among the first 2d — l lines of ^ d +i 
and then, one has no choice for the remaining: line and the columns: ( 2 ^) = ( 2d ~ 1 ) 
choices. 



To pick <^+i,z/Cd+i,2d+i for a given 2 < I < 2d, one needs to pick d — 1 lines among 
the first 2d — 1 lines of fi d +i and one cannot pick the (k — l)-th line. Then one has no 
choice for the remaining line and the columns: ( 2 d Zi) choices. 



d+l,2d+l 



We factorize the term 1 / 
and get the claimed bound. 

4 ) Proof of Proposition [751 - According to Proposition [3, 
-1 



whose log-expectation cancels out with E 7r0i ^ d log 



Cod 



log 2 



log p + E log 



Cod 



+ 



max (Elog |det ; Elog |det ; jty 1 ) + 7(^ 2 )) 



where 



Therefore, 



-1 
bg~2 



._a_ o 

Ota; 



and 



max (2Elog |ai| ; 2Elog |ax| + 2 log a ; 



2Elog|ai| +loga + 7(^ 1 ) +7(^ 2 )). 



Since a < 1, log a < 0, therefore 

-1 



log 2 



max 2Elog|ai| ; 2Elog|ai| +7 h/>M + 7(^ 2 ) 



(42) 
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where 

I _a o 

In order to finish the proof, we will construct of family of functions (/ p ) p gN* from (0, 1] x [0, 1] 
to R such that for all p £ N*, / p (a, /3) is non-decreasing in a and in (3 and such that for all 
PEN* and for all (a,/3) £ (0, 1] x [0, 1], 

7 (?) +tW> 2 ) </„(«, £)• 
We define T> in the following way: 

V 4 |J {(«, 0) £ (0, 1] x [0, 1] ; /„(a, /3) < 0}. 



Since for all p £ N*, / p (a, /3) is non-decreasing in a and in /?, we get that for all (x, y) £ Z>, 
(0,x) x [0,y) C V. Moreover, by (02J, if (a, (3) £ £>, then $l6j) is verified. 
Fix p £ N*. First note that by fl22|) . 

t(^) +7(^ 2 ) <-(Elog||^---^ 1 ||+Elog||^ 2 ---V 2 ||). 

Recall that we use the Frobenius norm on matrices. Denote </> 1 (a,/3) = 4 J l'''4 J l an d 
4> 2 (a,P) = ifjp • ■ -t^f. Note that the coefficients of and (j) 2 (a,f3) are polynomials in 

a and /3. The function 1/p (Elog ||</> 1 (o:, (3)\\ + Elog ||0 2 (a,/3)||) would be a good candidate 
for f p but it is not non- decreasing in a and (3, therefore, we have to modify it slightly. 
Consider P a polynomial in a and j3, 

n 

P(a,/3) = ^ ,,'>'•>'' • 

»>J=1 

Define the polynomial |P| in the following way 

« 

|P|(a,/5) = ^|^|a^. 

By the triangle inequality, for all (a, (3) £ (0,1] x [0,1], \P(a,/3)\ < \P\(a,(3). Moreover, 
\P\ (a, (3) is non decreasing for (a, (3) £ (0, 1] x [0, 1]. 

Define the matrices (a, /3) and |0 2 | (a, (3) in the following way. 
For i,j,k = 1,2, set \(f) k \^ = Then, 

I^V,/?)!! < || (a,/3)\\ and ||0 2 (a,/3)|| < || |0 2 | (a, /3)\\ . 
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Moreover || and || |0 2 | are non decreasing for (a, (3) G (0, 1] x [0, 1]. Thus, 

we conclude the proof by defining 

f p = l - (Elog || \<p l \ {a,(3)\\ + Elog || |0 2 | {a,(3)\\) . 



Remark 43 Note that if we define 



P\bi\ a \ I P\bj\ oK| 

|tti| I / L/.2I _ I \ci\ \&i 

\Ci\ 



and \ip i 






\ \<H\ 

then 

f p = 1/p (Elog I |^| • • • |$| J) + Elog || |^| • • • |^ 2 | ||) . 
We use that fact in the numerical computation of the functions f p . 

C. Exterior product 

In this section we give the material on exterior products. We provide only the properties 
relevant to the article, see [22, Chapter XVI. 6- 7] and [13, Chapter A. III. 5] for more details. 

Proposition 44 For < k < p, the exterior product of k vectors in C p , v±, . . . , is denoted 
by V\ A ■ • ■ A Vk- Is is a vector of the exterior product of degree k of C p that we denote by 
is a C-vector space of dimension ( fc ) . 
The exterior product v%, . . . ,Vk is a multi-linear (i.e. linear in every V{, 1 < i < k) and 
anti-symmetric (i.e. i^m A •••f CT (ifc) — e {°~) f or ° permutation of {l,...,k} and e{a) its 
signature) function. 

Ifei, . . . , e p is a basis ofC p , then (e^ A • • • Aejjl < %\ < ■ ■ ■ < < p) is a basis of f\ k C p . 

The later is called the canonical basis of f\ k C p ifei, ■ ■ • , &p is the canonical basis ofC p . 

If M is a matrix of size p x q, the exterior product of M that we denote by f\ k M is a map 

from f\ k C q to A fc C p such that 

k 

A M (vi A • • • A v k ) = Mvi A ■ ■ • A Mv k . 
Finally, for two matrices M and N, /\ k (MN) = f\ h (M) f\ h (N). 

Proposition 45 If X is a square matrix of size p, then 

v 

/\X = detX. 
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Moreover 



p 



detf\X = (detXf 



D. Capacity of the non-fading channels 

In this Section, we give expressions of the limiting sum-rate per-cell capacity for the Soft- 
Handoff model and the Wyner model (both symmetric and asymmetric) for the non-fading 
channels. Those expressions are consequences of results on Toeplitz matrices [23] . See [3] for 
an example of derivation. 

1 ) The Soft-Handoff model: We assume that d — 1, and for % G N*, Cm+i = E [0, 1] and 
Q t i = 1. Then, the limiting per-cell sum-rate capacity is 



2) The Wyner model: 

The symmetric setting: We assume that d = 2, and for % G N*, Cm+2 = Cm = 6 [0, 1] 
and Cm+i = 1- Then, the limiting per-cell sum-rate capacity is 

Cap(p) = [ log (l + Kp (1 + 2a cos(2vr/)) 2 ) df. 
Jo 

The asymmetric setting: We assume that d — 2, and for % G N*, Cv+2 — Ci,i+i — a ^ 
[0, 1] and Ci,i = 1- Then, the limiting per-cell sum-rate capacity is 
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Fig. 3. Comparison of the bounds of Corollary 1 1 1 1 and Proposition 1121 for A = 0.1, 1, in function of d. 
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Fig. 4. Comparison of the different bounds for A = 0.1, 1 and d = 3. 
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Fig. 5. Comparison of the different bounds for A = 0.1, 1 and K = 1,2, 4, 10. The values for the non random channel 
with A = 1 (resp. A = 0.1) and K = 1, 2, 4, 10 are respectively 1.06, 1.47, 1.95, 2.66 (resp. 2, 66, 3.25, 3.87, 4.72). 
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Fig. 7. Influence of the correlation on the capacity in function of the SNR and K . Bounds of Proposition [13] 
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Fig. 8. Influence of the correlation on the capacity in function of a. Bounds of Proposition 1 131 
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Fig. 9. Comparison of the symmetric and asymmetric channels with and without Rayleigh fading for A = 1, 10 4 in 
function of a. Bounds of Proposition 1131 
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Fig. 10. Comparison of the symmetric and asymmetric channels with Rayleigh fading for A = 10 4 and a = 0.5. 
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Fig. 11. Comparison of the symmetric and asymmetric channels with uniform fading for A = 1, 10 4 and a — 0.5 in 
function of e. Bounds of Proposition 1131 
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